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Background: Heavy-ion fusion reactions involving heavy nuclei at energies around the Coulomb barrier exhibit 
fusion hindrance, where the probability of compound nucleus formation is strongly hindered compared with that 
in light- and medium-mass systems. The origin of this fusion hindrance has not been well understood from a 
microscopic point of view. 

Purpose: Analyze the fusion dynamics in heavy systems by a microscopic reaction model and understand the 
origin of the fusion hindrance. 

Method: We employ the time-dependent Hartree-Fock (TDHF) theory as a microscopic reaction model. We ex¬ 
tract nucleus-nucleus potential and energy dissipation by the method combining TDHF dynamics of the entrance 
channel of fusion reactions with one-dimensional Newton equation including a dissipation term. Then, we analyze 
the origin of the fusion hindrance using the properties of the extracted potential and energy dissipation. 

Results: We obtain finite extra-push energies for heavy systems from TDHF simulations, which agree with 
experimental observations. Extracted nucleus-nucleus potentials show monotonic increase as the relative distance 
of two nuclei decreases, which induces the disappearance of an ordinary barrier structure of the nucleus-nucleus 
potential. This is different from those in light- and medium-mass systems and from density-constraint TDHF 
calculations. Extracted friction coefficients show sizable energy dependence and universal value of their magnitude, 
which are rather similar to those in light- and medium-mass systems. Using these properties, we analyze the origin 
of the fusion hindrance and Hnd that contribution of the increase in potential to the extra-push energy is larger 
than that of the accumulated dissipation energy in most systems studied in this article. 

Conclusions: We find that the nucleus-nucleus potentials extracted in heavy systems show a specihc property, 
which is not observed in light- and medium-mass systems. By the analysis of the origin of the fusion hindrance, 
we conclude that, as the system becomes heavier, the dynamical increase in nucleus-nucleus potential at small 
relative distances plays a more important role than the dissipation during the fusion reaction for understanding 
the origin of the fusion hindrance. 

PACS numbers: 21.60.Jz, 25.70.Jj 


I. INTRODUCTION 

Heavy-ion fusion reactions at energies around the 
Coulomb barrier have been studied for a long time and 
give interesting phenomena such as an enhancement of 
fusion cross sections at energies below the Coulomb bar¬ 
rier. For a better understanding of fusion dynamics, it 
is crucial to consider the couplings between the collec¬ 
tive motions of colliding nuclei and their internal exci¬ 
tations. To deal with such couplings, coupled-channels 
models have been developed and succeeded in reproduc¬ 
ing such an enhancement of fusion cross sections at sub¬ 
barrier energies In those models, unless a collid¬ 

ing system is too heavy or too light, it has been gen¬ 
erally considered that fusion takes place once colliding 
nuclei overcome the Coulomb barrier because of strong 
absorption inside the Coulomb barrier. Then, instead 
of using an imaginary potential inside the Coulomb bar¬ 
rier with the regular boundary condition at the origin, 
an incoming wave boundary condition Q has often been 
employed. The approaching phase leading to overcoming 
the Coulomb barrier to make a contact of projectile and 
target is considered as the capture process. 


The property of low-energy fusion dynamics drastically 
changes when the product of the charge numbers of the 
projectile and target nuclei, denoted by ZpZt, becomes 
greater than 1600. In such heavy systems, it has been 
known that the fusion probabiliW is extremely hindered 
around the Coulomb barrier [6|-l9|, which is called fusion 
hindranc4il- One of the well known experimental observa¬ 
tions for exhibiting this fusion hindrance is a direct com¬ 
parison between the -f ^®°Hf [ZpZp = 1296) and 
96Zr-f i24sn {ZpZp = 2000) systems leadi^ to ^^OTh 
reported in Ref. Q. The authors of Ref. Q measured 
evaporation-residue cross sections of these systems to de¬ 
duce their fusion probabilities for central collisions, and 
then compared them at the Coulomb barrier energy of 
both reactions estimated by the Bass model [a. They 
reported that the fusion probability of the ®®Zr-|-^^^Sn 
system is hindered by two orders of magnitude from that 
of the "^°Ar -f ^®°Hf system and from a simple barrier 
penetration model and that the apparent fusion barrier 
for the ®®Zr -f ^^^Sn system is shifted by about 26.7 MeV 
from the Bass barrier. Before this observation, Swiate- 
cki and coworkers had predicted that an additional ki- 


4 Notice that this phenomenon is different from fusion hindrance 
observed at deep subbarrier energies in light- and medium-mass 
systems [I3S1. 
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netic energy to the Coulomb barrier, called extra-push 
energy, is needed for fusion to take place when the so- 
called fissility parameter of the system (its definition 
given in Eqs. m and (HU) in Sec. IIII All is larger than 
its threshold value, which approximately corresponds to 
ZpZr > 1600 — 1800 [l3 - [r^ . They considered that, after 
the contact of colliding nuclei, large Coulomb force might 
make the colliding system reseparate under the appear¬ 
ance of an additional potential barrier inside the contact 
configuration. This indicates the necessity of treating the 
second step towards fusion, i.e., after the capture process 
as the first step, the process of shape evolution from the 
contact configuration of colliding nuclei to the configura¬ 
tion of compound nucleus formation. Since the projectile 
and target densities are significantly overlapped in this 
formation process, strong energy dissipation is expected 
to arise from the collective motions to the internal excita¬ 
tions. Therefore, the potential energy landscape and the 
mechanism of energy dissipation inside the Coulomb bar¬ 
rier are important quantities to understand the property 
of the fusion mechanism in heavy systems. 

For the analysis of the competition between compound 
nucleus formation and quasifission in such heavy systems, 
especially in the context of the synthesis of superheavy 
elements [isl - f^ . dynamical models such as a diffusion 
model based on a master equation [22l - [23| and a macro¬ 
scopic fluctuation-dissipation model based on a Langevin 
equation psl - f^ have been applied. In the latter model, 
relevant collective degrees of freedom such as the rela¬ 
tive distance and mass partition of colliding nuclei and 
their deformations are chosen and non-collective degrees 
of freedom are treated as heat bath. The inertia parame¬ 
ter, nucleus-nucleus potential, and friction coefficient for 
the relevant collective degrees of freedom, appearing as 
transport coefhcients in the equation, are important in¬ 
gredients for such calculations and are mainly determined 
by macroscopic models. 

The main motivation of the present article is to an¬ 
alyze the fusion hindrance in heavy systems by micro¬ 
scopic reaction models. In this work, we employ the 
time-dependent Hartree-Fock (TDHF) theory based on 
the Skyrme energy density functionals [29l - [^ as a micro¬ 
scopic reaction model. Energy density functional theory 
is selfconsistent, fully microscopic and provides a unique 
tool for describing both the static and dynamical proper¬ 
ties of nuclei across the whole nuclear chart in a unified 
framework [s^. Since the first applications of TDHF 
simulations to nuclear fusion reactions were per¬ 

formed about 40 years ago, much progress of TDHF sim¬ 
ulations has been achieved [s^ . Recent three dimen¬ 
sional TDHF calculations are able to include full Skyrme 
functionals consistent with nuclear structure calculations 
and have been ^plied to fusion reactions jiol - H^ , charge 
equilibration [^, transfer reactions as well as 

nuclear responses Furthermore, TDHF simula¬ 

tions to reactions with heavy nuclei have been performed 
[ssl - f^ . In Refs. [H, extra-push energies in 
several heavy systems have been estimated and the fu¬ 


sion hindrance resulting in finite extra-push energies has 
been realized, which is consistent with experimental ob¬ 
servations. In Refs. [13,[Ml, using the so-called density- 
constrained TDHF (DC-TDHF) method [i^, nucleus- 
nucleus potential as the minimized energy and excitation 
energy of a compound system are estimated along a path 
obtained from TDHF evolutions. Quasifission process is 
extensively studied in Refs. [5^,|^, 651, which have clari¬ 


fied the capability of TDHF simulations for the dynamics 
of quasifission. 

In the present article, we analyze the fusion hindrance 
in detail from the point of view of energy dissipation 
and dynamical nucleus-nucleus potential in the entrance 
channel of fusion reactions in heavy systems. For this 
purpose, we employ a method to simultaneously extract 
dynamical nucleus-nucleus potential and energy dissipa¬ 
tion for relevant collective coordinates from TDHF tra¬ 
jectories, which we proposed in Refs. [13, and called 
dissipative-dynamics TDHF (DD-TDHF) method. The 
DD-TDHF method relies on the assumption that a com¬ 
plicated microscopic mean-field dynamics in the entrance 
channel of fusion reactions can be reduced to a macro¬ 
scopic equation of motion with a dissipation term. Then, 
we apply this method to fusion reactions of heavy, nearly 
symmetric systems and analyze the property of nucleus- 
nucleus potential and energy dissipation extracted di¬ 
rectly from TDHF simulations. Using the properties of 
the extracted potential and dissipation, we perform a sys¬ 
tematic study of the analysis of the fusion hindrance in 
heavy systems. 

This article is organized as follows. In Sec. [Hi we ex¬ 
plain our method of how to analyze the fusion hindrance 
by TDHF simulations including how to extract nucleus- 
nucleus potential and energy dissipation from TDHF dy¬ 
namics. In Sec. imi we present the results of our system¬ 
atic calculations for extra-push energies, nucleus-nucleus 
potentials, energy dissipation in heavy systems, and the 
analysis of the fusion hindrance. We give the conclusions 
of the present article in Sec. lYl We explain some details 
of our method in the Appendix. 


II. METHOD 


A. Some remarks on TDHF 

In the TDHF theory, quantum many-body dynamics 
is replaced with one-body density dynamics, leading to 
the TDHF equation 

= [h[p], p], (1) 

where p and h[p] denote the one-body density and self- 
consistent single-particle Hamiltonian as a functional of 
one-body density, respectively. The selfconsistent single¬ 
particle Hamiltonian can be obtained from the energy 
density functional, denoted by S[p\, through a variation 
h[p\ = d£[p\/dp. 
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Since the TDHF theory describes collective motions in 
a semiclassical way, one can define semiclassical trajec¬ 
tories corresponding to the time evolutions of relevant 
collective coordinates from one-body density. In this ar¬ 
ticle, we concentrate on low energy nuclear dynamics, 
where the energy of the system is low enough to vali¬ 
date the use of the mean-field dynamics and is above the 
Coulomb barrier. 


B. DD-TDHF method 

In this subsection, we briefl y r ecall the DD-TDHF 
method, proposed in Refs. [3, The method is de¬ 
scribed as follows. First, we properly define the collective 
coordinates to describe fusion reactions. In this article, 
we concentrate on central collisions of fusion dynamics 
and employ the relative distance R between projectile 
and target nuclei as a collective coordinate. An extension 
to using several collective coordinates is possible. Then, 
we construct two-body kinematics of projectile and tar¬ 
get nuclei. In TDHF simulations, we define at each time 
a separation plane to separate the total density to the 
densities of two subsystems. The plane, called neck, is 
perpendicular to the collision axis, which is always a;-axis 
in the case for central collisions. We will explain how to 
determine the neck position in Appendix The densi¬ 
ties of the two subsystems are computed in the coordinate 
space representation as 

Pi, 2 (t', t) = p{r, t) e[±{x - a;o(t))], (2) 

where 0{x) is the step function and the neck position 
is at x = xo{t) depending on time. Then, we com¬ 
pute the center-of-mass coordinate Ri and associated 
momentum Pi as well as the mass number Ai for each 
subsystem i = 1,2 from the density of each subsys¬ 
tem, A^{t) = Tr(/5i(t)), R,{t) = TT{xpi{t))/Ai, and 
Pi{t) = Tr{pxPi{t)). The masses of the subsystems are 
directly computed as 


rrii = Pi/Ri, (3) 

for i = 1,2. From these quantities, we construct a TDHF 
trajectory for the relative distance R, associated momen¬ 
tum P, and reduced mass p at each time by 


R — Ri — i? 2 ) 

p _ m 2 Pi - TO1P2 
mi -I- TO2 
mim2 

M ^-• 

mi + m2 


( 4 ) 

( 5 ) 

( 6 ) 


We assume that the time evolutions of R and P obtained 
from the TDHF evolution of central collisions follow a 
one-dimensional macroscopic equation of motion with a 
dissipation term that is assumed to depend on the veloc¬ 


ity it: 


dR _ P 
dt p ’ 

dP _ dV d /P 2 \ dR 
dt dR dR V 2/r y ^ dt ^ 


( 7 ) 

( 8 ) 


where V and 7 denote the nucleus-nucleus potential 
and friction coefficient, respectively. This assumption 
is based on the hypothesis that a complicated micro¬ 
scopic dynamics obtained from the TDHF evolution can 
be reduced to a simple one-dimensional Newton equation. 
The aim is to obtain transport coefficients in the equa¬ 
tion from TDHF dynamics. We have already obtained 
the reduced mass by Eq. ®. To obtain two unknown 
quantities dV/dR and 7 as a function of R, we need two 
equations. If we assume that the transport coefficients, p, 
dV/dR, and 7 , at each R do not change by a slight change 
of energy AE, we can construct a system of two equa¬ 
tions for Eq. (jS]) corresponding to two TDHF trajectories 
at two slightly different energies Ej and Eu = Ej + AE. 
By solving the system of two equations at R, we obtain 
7 (i?) and dV/dR as 


i{R) = 

dV{R) 
dR ~ 


Pi — Pii Ri + Rii dp 
~Ri-Rii ^ 2 dR’ 

itjPjj — itjjPj RjRii dp 

Ri - Rii ~dR’ 


(9) 

( 10 ) 


where the subscripts I and II correspond to the trajecto¬ 
ries at energies Ei and Eu, respectively, and we perform 
all the time derivatives at R/ = R 77 = R. The nucleus- 
nucleus potential V is obtained from Eq. m by integra¬ 
tion with the asymptotic form of the Coulomb potential 
at sufficiently large distances. We employ Rmax = 20 fm 
in this article. Therefore, we do not need any normal¬ 
izations for the potential. In Refs. [13,14^, we applied 
the DD-TDHF method to fusion reactions in light- and 
medium-mass systems and gave the following results: (i) 
At energies near the Coulomb barrier, dynamical effects 
significantly decrease the barrier height of the extracted 
nucleus-nucleus potentials because of reorganization of 
the density of colliding nuclei, resulting in energy depen¬ 
dence of potential, (ii) The extracted friction coefficients 
show a universal property similar to that in one-body 
dissipation models, and significant energy dependence re¬ 
lated to the picture of the window formula of one-body 
dissipation mechanism. From those results, we concluded 
that the DD-TDHF method is useful tool to investigate 
the property of nucleus-nucleus potential and one-body 
energy dissipation in low-energy nuclear reactions. 

When we apply the DD-TDHF method to fusion re¬ 
actions in heavy systems, we find that, as the overlap 
of the densities of colliding nuclei becomes substantially 
large at small relative distances in TDHF simulations, 
constructing two-body kinematics of colliding nuclei with 
a proper determination of the neck position becomes dif¬ 
ficult. Also, we face in some systems the following situ- 
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ations: (i) the relative velocity or the velocity of a frag¬ 
ment becomes 0 , (ii) the mass of a fragment becomes 
negative. We consider that these situations also make it 
impossible to construct two-body kinematics. Therefore, 
we stop performing the extraction once we face these sit¬ 
uations or once we can not define the neck position. In 
the following, we denote as i?niin the relative distance at 
which we stop the extraction. 


C. Nucleus—nucleus potential from the frozen 
density approximation 

As a reference for the nucleus-nucleus potential based 
on a framework of the energy density functional, we em¬ 
ploy the so-called frozen density approximation [^, . 

In this approximation, the energy of the total system 
at a given relative distance R between the centers-of- 
mass of the projectile and target densities, denoted by 
£[pp+t\{R), is obtained with keeping their densities 
frozen to their respective ground-state densities. Denot¬ 
ing the ground-state densities of the projectile and target 
nuclei as pp and pp, respectively, we obtain the frozen 
density potential as 

= SIpp+tKR) - £[pp] - £Ipt], (11) 

where £[p\ is calculated from the same Skyrme energy 
density functionals as those used in our TDHF simula¬ 
tions. In the frozen density potential, effects of possible 
rearrangement of internal structure during collision are 
not considered. Note that we neglect the Pauli effect from 
the overlap of the projectile and target densities in esti¬ 
mating £[pppt\(R)- Namely, we approximate pp+p as a 
simple sum of pp+pp. This will be crucial when the over¬ 
lap becomes large at small relative distances, leading to 
overestimation of the attractive potential energy. There¬ 
fore, we use the frozen density potential as a reference 
until a relative distance around the Coulomb barrier ra¬ 
dius, where the overlap of projectile and target densities 
is expected to be small. 


D. Definition of fusion in heavy systems 

Semiclassical nature of TDHF dynamics enables us to 
define a threshold energy for classification of nuclear re¬ 
actions. In the case of central collisions at low energies, 
we can define the fusion threshold energy above which 
fusion, i.e., the formation of a compound nucleus, takes 
place. Following a manner similar to those in previous 
works m, 113, we define fusion as a reaction where a 
colliding system keeps a compact shape for a sufficiently 
long time (~ 2000 fm/c) after contact. From a sufficiently 
low energy, we perform TDHF calculations systemati¬ 
cally by increasing progressively the energy with a step 
of 0.5 MeV until fusion is reached. 

We have found that there is an ambiguity to define a 
compact shape leading to the formation of the compound 


nucleus in some systems. We will come to this point in 
Appendix [HI and will show that this ambiguity does not 
affect the conclusion of this article. 


E. Numerical setups 


In this work, we perform only central collisions. In 
our practical computations of TDHF evolutions, we use 
the TDHF3D code [13 developed by P. Bonche and 
coworkers with the SLy4d parametrization [43 of the 
Skyrme energy density functional, which is slightly mod¬ 
ified from the SLy4 parametrization Hi, [73 by neglect¬ 
ing the center-of-mass corrections in the fitting protocol 
of determining the parameter set. As the initial condi¬ 
tions, we solve static Hartree-Fock equations by using the 
computer code ev 8 74| with the imaginary time method. 
Then, we place the centers of the projectile and target 
densities on the a;-axis. The initial distance is set to be 
Rq = 28.8 fm for all the systems considered in this article 
except for the "‘°Ca-|-"^°Ca system {Rq = 22.4fm). We 
assume that the colliding nuclei follow the Rutherford 
trajectory before the initial distance. The step sizes in 
the coordinate space and time are 0.8 fm and 0.45 fm/c, 
respectively, for both static and dynamical computations. 
The numbers of the mesh points in our three dimensional 
numerical box is 64 x 28 x 14 along the x, y, and z-axis, re¬ 
spectively, with a symmetry plane at z = 0. The pairing 
correlation is neglected in both static and dynamical cal¬ 
culations, though this has been included in some of recent 
applications to nuclear responses and reactions in either 
time-dependent Hartree-Fock-Bogoliubov [sE HE HE or 
an approximated way [IE, [ 


III. RESULTS 


To analyze the fusion hindrance, we perform TDHF 
simulations for central collisions in heavy systems 
90 zj.-hW 92 . 94 . 962 i., iooMo-H ®2’i™Mo,i°4Ru,“°Pd, and 
+ 90.92.94.96Z]., where the extra-push energies were 
deduced from experimental evaporation-residue cross sec¬ 
tions [ 1 , 13 , and 96Zr-Hi32gjj,i36xe, 70zn + 208 pb^ 

We obtain the ground-state shapes of ®^Zr, ®^Mo, 
ioomo, and ^^"^Sn to be oblate with the deformation pa¬ 
rameters p = 0.03, 0.05, 0.23, and O.ll, respectively. The 
ground-state shapes of ®^Zr and ^^®Xe are slightly pro¬ 
late deformation with /3 = 0.04 and 0.05, respectively. 
The ground-states of ^°^Ru and ^^°Pd are triaxial with 
P = 0.31, 7 = 21.0° and P = 0.31, 7 = 21.6°, respec¬ 
tively. The other nuclei are found to be spherical. Note 
that the ground-states of ®^’®^Zr, ®^Mo, and ^^®Xe 

would be spherical if the pairing correlation were included 
in the calculations. For collisions with deformed nuclei, 
we place the nucleus in the box by setting the symmetry 
axis of the axially deformed nucleus or the longest axis 
of the triaxial nucleus parallel to the z-axis. 
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FIG. 1. (Color online) Time evolution of mass quadrupole 
moment Q 20 of the total system at different center-of-mass 
energies for the ®°Zr + ^^“^Sn system. 
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FIG. 2. (Color online) Time evolution of density profiles 
p{x,y,z = 0,t) at Ec.m. = 233.5 MeV for the ®°Zr + '^^'‘Sn 
system. 


A. Extra-push energy 

In our TDHF simulations in heavy systems, we indeed 
observe that, after contact of colliding nuclei, a collid¬ 
ing system keeps a configuration of substantial overlap of 
projectile and target densities for a long time, and then 
reseparates to two fragments. We regard this process as 
quasifission. This has been observed in previous works 
[ 3 ^ and we are not aware of such a process in 

lighter systems. As an example of quasifission realized in 
our TDHF simulations, we show in Fig. [T] the time evo¬ 
lution of mass quadrupole moment of the ®*^Zr-|-^^^Sn 
system defined as 

Q 2 o{t) = J (fr{2x'^-y^ - z^)p{x,y,z,t), (12) 

where p{x, y, z, t) denotes the density of the total system 
obtained from TDHF calculations at different center-of- 
mass energies Ac.m. In this system, we can identify that 
the case at E^.m. = 233.5 MeV leads to quasifission af¬ 
ter a long contact time ^ 1800 fm/c. This energy is well 
above the barrier height obtained from the frozen den¬ 
sity potential, 217.1 MeV. In Fig. [2J we show the time 


evolution of density profile for the quasifission process 
at Ec.m. = 233.5 MeV. Large overlap of the densities 
from t = 450 fm/c to t = 2250 fm/c and reseparation 
at t = 3015 fm/c can be seen. 

We define the extra-push energy Aextra estimated from 
TDHF simulations as 


Aextra — Athres VpD 5 


(13) 


where Flthres and Vfd denote the fusion threshold energy 
and the potential barrier height obtained from the frozen 
density approximation, respectively. Figure [3] shows the 
extra-push energies obtained from TDHF simulations de¬ 
noted by the red filled circles for various heavy systems. 
We also show by the blue filled triangles experimental 
extra-push energies taken from Ref. @ , which are defined 
as the difference between the apparent barrier deduced 
from experimental data of evaporation-residue cross sec¬ 
tions and the barrier height obtained from the Bass model 
Our TDHF simulations reasonably reproduce the 
experimental extra-push energies, except for heavier sys¬ 
tems. We also reproduce the trend that the extra-push 
energy becomes large as the charge product of the sys¬ 
tem becomes large. Figure 3] compares our result with 
the extra-push model [lil - ll^ , plotted by the dot-dashed 
line. In this figure, we plot extra-push energy as a func¬ 
tion of parameter Xmea.n defined as 

1 2 

^mean — ^^efF “b (I^) 

with 


4ZpZt / 

( 4/3 + 4 / 3 ) / 

xns=(Z^/A)/(Z^/Al^.^^. 


(15) 

(16) 


Here, (Z^/A)crit = 50.883(1 — 1.7826/^) is used, where 
I = (A — 2Z)IA is the neutron excess of the compound 
nucleus with A = Ap + Ap, Z = Zp + Zp being the 
mass and charge numbers of the compound nucleus, re¬ 
spectively. The parameter Xmean was originally proposed 
in Ref. [13, with which the extra-push energy better 
scales. Around a^mean = 0.74, our result is better than 
the extra-push model, while the extra-push model is bet¬ 
ter at a^mean ^ 0.77. 

In some systems, we find that our TDHF simulations 
underestimate the values of the experimental extra-push 
energy. We underestimate the extra-push energy by 
about 12 MeV in the ®®Zr-|-^^^Sn system. A possible 
interpretation of this underestimation is as follows. As 
we have explained in Sec. IHD[ we define fusion by TDHF 
as a reaction where a colliding system keeps a compact 
shape for a sufficiently long time. As one can see in 
Fig. [2 such a reaction never reaches to the shape of 
resulting compound nucleus in TDHF simulations be¬ 
cause of lack of dissipation beyond the mean-held level, 
e.g., two-body dissipation arising from direct nucleon- 
nucleon collisions [72 h 74| |. Therefore, the cases dehned 
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FIG. 3. (Color onlinelExtra-push energies extracted from TDHF simulations (red filled circles) and from experimental data 
(blue filled triangles) [y]. In the a;-axis, the nuclide of the target is shown, while the nuclide of the projectile is shown at the 
top. 



^mean 

FIG. 4. (Color online) Same as in Fig.[3l but as a function of 
parameter Xmean defined by Eq. (fT4t . 


as fusion from TDHF simulations might lead to resepara¬ 
tion (quasifission) if dissipation mechanisms beyond the 
mean-field level were included in the calculations. That 
is, fusion threshold energies in some systems might be 
underestimated within the mean-field level. Other possi¬ 
ble interpretation of this underestimation is the neglect 
of pairing correlations in both static and dynamical cal¬ 
culations. This might induce the discrepancy between 
TDHF and experimental fusion threshold energies. 

Although we have such a limitation, our TDHF sim¬ 
ulations qualitatively reproduce the overall property of 
the dependence of observed extra-push energies on sys¬ 
tem size. This fact validates our following analysis based 
on TDHF simulations. In the following, we will give de¬ 
tailed analyses of nucleus-nucleus potential and energy 
dissipation obtained from the DD-TDHF method and of 
the origin of the fusion hindrance in heavy systems. 



R [fm] 

FIG. 5. (Color online) Nucleus-nucleus potential as a function 
of relative distance R for the ^^Zr + ^^'^Sn system. The red 
solid, green dot-dashed, blue dashed, and brown dotted lines 
denote the potentials extracted from the TDHF trajectories 
at -Ec.m. = 300, 250, 228, and 220 MeV, respectively. The 
filled circles are the potential obtained from the frozen density 
approximation. 


B. Nucleus—nucleus potential in heavy systems 

1. Property of extracted potential 

Figure shows the nucleus-nucleus potentials as a 
function of relative distance R obtained from the DD- 
TDHF method for the ®®Zr-|-^^“^Sn system at i?c.m. = 
300, 250, 228, and 220 MeV used in TDHF. The filled 
circles are the potential obtained by the frozen density 
approximation. Note that, when R < 11.2 fm, the over¬ 
lap density of the two nuclei at the neck in the frozen den¬ 
sity approximation becomes close to the normal density 
0.16fm“^, indicating the violation of the frozen density 
approximation at i? < 11.2 fm. The charge product of 
this system is 2000, and we observe the fusion hindrance 
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FIG. 6. (Color online) Same as in Fig. [5l but for the 
'^°Ca + '^°Ca system. 

in this system as we have shown in Fig. |3l The fusion 
threshold energy for this system is E^.m. = 228 MeV. At 
Ec.m. = 220 MeV, which is higher than the barrier height 
of the frozen density potential, the system reseparates af¬ 
ter slight contact. Note that, for the cases at E^.m. = 250, 
228, and 220 MeV, we stop the extraction of the potential 
at i?min = 10.1, 11.4, and 12.7fm, respectively. The ex¬ 
tracted potentials agree with the frozen density potential 
until R 13.5 fm and start to deviate from the frozen 
density one as R decreases. At R < 13.5 fm, the ex¬ 
tracted potentials do not have an ordinary barrier struc¬ 
ture, though the frozen density potential has an ordinary 
barrier at i? ~ 12.8 fm. The extracted potentials mono- 
tonically increase after deviating from the frozen den¬ 
sity potential. At i? > 11.4 fm, the extracted potentials 
agree with each other, while a deviation is seen between 
the potentials extracted at Ec.m. = 300 and 250 MeV 
at i? < 11.4 fm. To investigate these findings in detail, 
we compare them with the potentials extracted in light 
systems. Figure H] show the potentials extracted in the 
^°Ca-|-^°Ca system. From the comparison, we find two 
main differences between the potentials extracted in light 
and heavy systems as follows: 

• Potentials for the '^°Ca-|-'^°Ca system and the 
frozen density potential for both systems have an 
ordinary barrier structure, while no barrier appears 
in the potentials extracted in the ®®Zr-|- ^^"‘Sn sys¬ 
tem. Monotonic increase of the extracted poten¬ 
tials is only observed in the ®®Zr -|- ^^'^Sn system. 

• In the ^^Ca. + ^^Ca. system, energy dependence of 
potential around the Coulomb barrier is significant, 
while for the ®®Zr-|-^^^Sn system, energy depen¬ 
dence of potential is less pronounced than in the 
40ca-p40ca system. 

The above two properties are indeed general in heavy 
systems. To see this, we show in Fig.Qextracted nucleus- 
nucleus potentials for the ^°°Mo -I- ^°'^Ru {ZpZt = 1848), 
loOMo-piiOpd (1932), 9°Zr-bi24sj^^ 962r + i32sn (2000), 
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FIG. 7. (Golor online) Nucleus-nucleus potentials for various 
heavy systems are displayed at the region around the barrier 
of the frozen density potential. 

96zr-pi36xe (2160), ™Zn-f208pb (2460) systems. In 
each system, except for the ™Zn + ^°®Pb system, we 
show the potentials extracted at the fusion threshold en¬ 
ergy by the blue dashed line and at higher energy by 
the red solid line. In the ™Zn-|-^°®Pb system, we do 
not observe fusion at energies between 265 and 400 MeV. 
In any systems, the potentials show monotonic increase 
as R decreases and energy dependence of potential is 
less pronounced. (At most, the difference is 0.6% in 
the ®°Zr-f ^^"^Sn system, while the difference is 3% in 
the "^°Ca-|-"^°Ca system.) By this systematic study of 
nucleus-nucleus potential, we conclude that these prop¬ 
erties of the extracted potentials are general in heavy 
systems. 

8. Discussion 

We have shown that the nucleus-nucleus potentials ex¬ 
tracted in heavy systems exhibit monotonic increase as 
the relative distance decreases and do not have an or¬ 
dinary barrier structure and this property is seen only 
in heavy systems. In the following, we will consider the 
reason why this property appears only in this case. 

Figure [5] shows the relative momentum scaled by the 
initial reduced mass number Aq = ApAt/{Ap -|- Ap) as 
a function of relative distance for the ®®Zr-|-^^^Sn (red 
solid line) and the ^°Ca-|-^°Ca (blue dashed line) sys¬ 
tems. The Coulomb barrier radius of the frozen density 
potential is indicated by arrow for these systems. Note 














FIG. 8. (Color online) Absolute value of the relative mo¬ 
mentum scaled by the initial reduced mass number Aq as a 
function of relative distance for the ®®Zr-|-^^^Sn at E^.m. = 
228 MeV (red solid line) and the ^°Ca-|-^°Ca at Ec.m. = 
60 MeV (blue dashed line) systems. The barrier radii of the 
frozen density potentials for these systems are indicated by 
arrow. 



FIG. 9. Schematic illustration of the appearance of a condi¬ 
tional saddle in heavy systems. See text for details. 


We further consider the meaning of the increase in po¬ 
tential at small relative distances. To do so, we introduce 
with Fig.j^an explanation of the appearance of the fusion 
hindrance in heavy systems from a simple geometrical 
feature of potential used in the extra-push model . 

Figure [S] shows the schematic illustration of an adiabatic 
potential energy landscape as a function of deformation 
of a compound nucleus in light and heavy systems. In 
both systems, the left region corresponds to the configu¬ 
ration of a compound nucleus formation, while the right 
corresponds to the quasifission. The potential barrier, 
denoted by “saddle” in the figure, is usually called condi¬ 
tional saddle and corresponds to the fission barrier under 
the condition that the mass partition of the system into 
two fragments is fixed to be the one at the contact con¬ 
figuration in the entrance channel. If the equilibration at 
the contact configuration is fast enough, in other words, 
the transition from two-body regime to one-body regime 
in the collision dynamics is fast enough, the system will 
feel this adiabatic potential once two nuclei touch. In 
light systems, the contact configuration is located inside 
the saddle on the potential landscape. Consequently, fu¬ 
sion automatically takes place once the two nuclei touch. 
The situation is changed for heavy systems because of 
large Coulomb energy. That is, the saddle is shifted to 
the inside of the contact configuration, and the system 
needs to overcome the saddle for fusion, which gives rise 
to the fusion hindrance in the extra-push model. In the 
DD-TDHF method, dynamical effects are automatically 
included in the extracted potential through TDHF dy¬ 
namics. Since the TDHF theory treats the reorganization 
of single-particle wave functions during collision in a self- 
consistent way, the extracted potential includes effects of 
smooth transition from two-body regime of collision dy¬ 
namics to one-body regime of shape evolution at small 
R. Therefore, we consider that the increase in potential 
at small R observed in heavy systems can be regarded as 
a reflection of the appearance of the conditional saddle 
inside the contact configuration in the schematic picture 
explained above. 


that the sign of the relative momentum is negative in the 
entrance channel of central collisions. In the "'^°Ca-|- ^°Ca 
system, the momentum decreases as R decreases until 
the Coulomb barrier, and then the momentum starts to 
increase after the system overcomes the Coulomb bar¬ 
rier. In the ®®Zr -I- ^^'^Sn system, however, the momen¬ 
tum monotonically decreases as R decreases. If the po¬ 
tential of this system were the frozen density one, the mo¬ 
mentum would increase after passing the barrier. There¬ 
fore, the monotonic decrease of the momentum with R 
decreasing induces the increase in extracted potentials 
at the inside of the frozen density barrier. Dissipation 
also can decrease the momentum. However, our anal¬ 
ysis shows that the extracted dissipation is not enough 
to explain the property of the time evolution of the mo¬ 
mentum in heavy systems and the increase in potential 
is necessary. 


3. Isotope dependence of nucleus-nucleus potential 

Next we will analyze the isotope dependence of 
nucleus-nucleus potentials in heavy systems. For this 
purpose, we employ the ®®Zr-|-^^^Sn and ®®Zr-|-^^^Sn 
systems by changing the neutron number of Sn. 

In Fig. fTUl a). we show the nucleus-nucleus potentials 
for the ®®Zr -|- ^^^Sn (red solid line) and ®®Zr -|- ^^^Sn (blue 
dashed line) systems. The potentials in both systems are 
extracted at E^.m. = 230 MeV. As a reference, the frozen 
density potentials are plotted by the red filled circles and 
blue filled triangles for the ®®Zr -|- ^^“^Sn and ®®Zr -|- ^^^Sn 
systems, respectively. In the DD-TDHF method and 
frozen density approximation, the potential energy for 
the ®®Zr-|-^^^Sn system is slightly larger than that for 
the ®®Zr-|-^^^Sn system. In the DD-TDHF case, the 
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FIG. 10. (Color online) (a) Comparison between the nucleus- 
nucleus potentials for the ^^Zr + ^^'^Sn (red solid line) and 
®®Zr + ^®^Sn (blue dashed line) systems. Both potentials are 
extracted at Sc.m. = 230 MeV. The red filled circles and 
blue filled triangles are the frozen density potentials for the 
^^Zr + ^^'^Sn and ®®Zr + ^®^Sn systems, respectively, (b) Av¬ 
erage number Atrans of transferred neutrons (red) and pro¬ 
tons (blue) from ^^“^Sn (solid line) and ^®^Sn (dashed line) 
to ®®Zr in the ®®Zr-f systems at Ac.m. = 230 MeV. 

(c) Neutron to proton ratio {N/Z) of the two subsystems in 
the ®®Zr -I- ^^"^Sn (solid line) systems and ®®Zr + ^^^Sn (dashed 
line) systems. The gray solid and dashed lines indicate the 
values of the equilibrated N/Z for the two systems. 


difference between the potentials is 0.4 MeV at the bar¬ 
rier radius of the frozen density potential, R = 12.8 fm, 
which is similar to the difference between the frozen 
density potential barriers of the two systems, 0.5 MeV. 
Then, the difference between the two extracted poten¬ 
tials evolves as the relative distance becomes small, and 
becomes 1.8 MeV at R = 11.1 fm. In order to further 
understand the role of the neutron excess in ^^^Sn on 
the potential, we show in Fig. lTUI b) the average numbers 
of transferred nucleons, Atrans, through the neck from 
i24Sn and to in the ^^Zr-h i24.i32gj^ systems 

at Ec.m. = 230 MeV. The positive and negative values 
correspond to transfer from to ®®Zr and trans¬ 
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FIG. 11. (Color online) Same as Fig. [in] but for the 
90,96 Zj.-1- i24gn systems obtained at the fusion threshold en¬ 
ergy ifthres = 234.5 MeV and ifthres = 228 MeV, respectively. 


fer from ®®Zr to respectively. We find that 

the number of transferred neutrons in the ®®Zr-|-^^^Sn 
system is larger than that in the ®®Zr-|- ^^“^Sn system as 
R decreases. We also find the larger number of trans¬ 
ferred protons from ®®Zr to ^^^Sn than that from ®®Zr to 
^^^Sn. Because of larger neutron to proton {N/Z) ratio 
of 132Sn {N/Z = 1.64) than that of i24Sn {N/Z = 1.48), 
larger nucleon transfer towards charge equilibration oc¬ 
curs during the collision in the ®®Zr -|- i32Sn system. This 
is illustrated in Fig. uni: c), where N/Z of the two sub¬ 
systems in the ®®Zr-|-i^^^Sn (solid line) and ®®Zr-|-i32Sn 
(dashed line) systems are plotted. As R decreases, larger 
change in N/Z towards the change equilibration in the 
®®Zr-|-i32Sn system can be seen. We consider that the 
net effects make the potential for the ®®Zr -|- 1^2 Sn system 
lower. 

Similar property is found in the extracted potentials 
in ®°’®®Zr -I- i^^Sn systems. Figure |TT[a) shows the ex¬ 
tracted potentials at the fusion threshold energy and 
the frozen density potentials for these systems. We find 
that the increase of the extracted potential from the 
frozen density potential is larger in ®®Zr-|-i^^Sn than 
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FIG. 12. (Color online) Comparison between the nucleus- 
nucleus potential obtained from our method (red solid line) 
and that from the DC-TDHF method (blue dashed line) 
taken from Ref. for the ®®Zr + ^^^Sn system at Ec.m. = 
230 MeV. As a reference, the frozen density potential of this 
system is plotted by the filled circles. 

in ®°Zr + ^^'^Sn. The difference is about 4 MeV in aver¬ 
age. From Fig. Illf b) showing the number of transferred 
nucleons and Fig. Illf c) showing N/Z ratio, larger nu¬ 
cleon transfer towards charge equilibration occurs in the 
®°Zr-|- ^^“^Sn system due to larger difference between the 
initial N/Z ratios of the projectile and target. 

From these observations, we conclude that, as nucleon 
transfer towards charge equilibration during the collision 
becomes smaller, increase in the extracted potential from 
the frozen density potential becomes larger. 


4- Comparison with the DC-TDHF method 

Figure [1^ compares nucleus-nucleus potentials ob¬ 
tained from our method (red solid line) and from the 
DC-TDHF method (blue dashed line) in Ref. [6l| for 
the ®®Zr-|- ^^'^Sn system at ifc.m. = 230 MeV, which is 2- 
MeV above the fusion threshold energy. As a reference, 
the frozen density potential is plotted by the filled cir¬ 
cles. The DC-TDHF potential has an ordinary barrier at 
R ~ 13.1 fm, which is slightly outside the barrier radius 
of the frozen density potential, R ~ 12.8 fm. The barrier 
height of the DC-TDHF potential is about 3.5-MeV lower 
than that of the frozen density potential. The difference 
between DC-TDHF potential and our potential becomes 
large as the relative distance becomes small. 

Since both methods solve the same TDHF equations, 
the time evolutions of the single-particle wave functions 
in the two methods should be same, except for numerical 
errors and slightly different Skyrme’s parametrizations 
(SLy4d for DD-TDHF and SLy4 for DC-TDHF). In the 
DC-TDHF method, the potential energy is given as the 
energy density functional of the total system with the 
minimization under the constraint that the total density 



R [fm] 

FIG. 13. (Color online) Same as in Fig. [12] but for the 
'‘°Ca-b'‘°Ca system at Sc.m. = 55 MeV. DC-TDHF potential 
is taken from Ref. (75l| . 

matches the TDHF density. Therefore, the obtained po¬ 
tential can be regarded as the lowest energy along the dy¬ 
namical path that includes all dynamical effects through 
the time evolution of the realistic TDHF density distribu¬ 
tions. In the DD-TDHF method, the potential energy is 
obtained under the assumption that complex TDHF dy¬ 
namics for central collisions is reduced to one-dimensional 
macroscopic equation of motion. Effects coming from 
other collective degrees of freedom than the degree of 
freedom of the relative motion are reduced to the trans¬ 
port coefficients of the reduced mass, nucleus-nucleus 
potential, and dissipation term in the one-dimensional 
Newton equation. Although both methods are based on 
the same TDHF dynamics, the way of interpretation of 
nucleus-nucleus potential is totally different from each 
other, giving different properties of the potential in heavy 
systems. Note that the reduced mass extracted from the 
DD-TDHF method includes dynamical effects, leading to 
the R dependence of mass. 

However, in light systems, we do not see a signifi¬ 
cant difference between DC-TDHF and DD-TDHF po¬ 
tentials. As an example, we show in Fig. |T2| the com¬ 
parison between DC-TDHF (blue dashed line) and our 
potential (red solid line) for the ^°Ca-|-^°Ca system at 
Ec.m. = 55 MeV, and find no significant difference be¬ 
tween the two potentials (only 0.3 MeV difference at most 
at the barrier). Therefore, the effects mentioned above 
are significant only in heavy systems. 

C. Dissipation 

The TDHF theory includes one-body dissipation 
through the selfconsistent mean field. As we explained 
in Sec. IHB[ we can simultaneously extract the nucleus- 
nucleus potential and the friction coefficient from the 
DD-TDHF method. In this subsection, we discuss the 
property of extracted friction coefficients in heavy sys¬ 
tems. 
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FIG. 14. (Color online) Reduced friction coefficient defined as 
7 (i?)/^(i?) as a function of R extracted from TDHF trajec¬ 
tories at Ec.m. = 300 MeV (red solid line), Ec.m. = 250 MeV 
(green dot-dashed line), i?c.m. = 228 MeV (blue dashed line), 
and Ec.m. = 220 MeV (brown dotted line) for the ®®Zr -|- 
system. 



FIG. 16. (Color online) Accumulated dissipation energy cal¬ 
culated by Eq. (O as a function of R obtained from TDHF 
trajectories at Ec.m. = 300, 250, 228, and 220 MeV for the 
®®Zr -I- ^^“^Sn system. A zoom around i? = 13 fm is inserted. 



FIG. 15. (Color online) Reduced friction coefficient 
^{R)/IJ.{R) as a function of R for selected systems. 


In Fig. [T31 we show the reduced friction coefficients 
defined as the friction coefficient divided by the reduced 
mass 7 (i?)//i(i?) as a function of relative distance in the 
®®Zr-|-^^“^Sn system extracted at the same energies as 
used in extracting potential in Fig. [5] The extracted fric¬ 
tion coefficient increases as the colliding nuclei approach 
at i? > 13 fm. And then, unlike the extracted poten¬ 
tial, we observe a clear energy dependence of the friction 
coefficient. As Ec.m. increases, the peak position of the 
extracted friction coefficient is shifted towards small R, 
which is close to i?min) the relative distance at which we 
stop the extraction, at each energy. (Note again that, for 
the cases at Ac.m. = 250, 228, and 220 MeV, i?min = 10.1, 
11.4, and 12.7fm, respectively.) The reason of the ap¬ 
pearance of the peak near Amin is that, as R becomes 
small, the relative velocity becomes small and the inter¬ 
nal structure of the colliding nuclei has much time to 


reorganize, giving a large value of the friction coefficient. 

Figure [15] shows the reduced friction coefficients for 
selected systems. For these systems, each friction co¬ 
efficient is extracted at each fusion threshold energy. 
The radial dependence of the friction coefficients is sim¬ 
ilar in these systems: As R decreases, the friction co¬ 
efficient increases, makes two humps at R ^ 12.6 and 
12 fm, and then has a sharp peak. Also, the values of 
their strengths at the humps are similar to each other, 
7 //r^ lxl0^^s“^, and their strengths at the sharp peak 
are of the same order of magnitude, 'j/fj, ^ 2—6x 10^^ s“^. 
For comparison, we also show the reduced friction coef¬ 
ficient for the "‘“Ca-I-"‘“Ca system at Ec.m. = 55 MeV by 
the open squares. The radial dependence of the friction 
coefficient is slightly different from those in heavy sys¬ 
tems: It has only a single peak. Its strength is however 
of the same order of magnitude as that at the peak in 
heavy systems. 

With the extracted friction coefficient, we can evaluate 
accumulated dissipation energy from the relative motion 
to internal excitations by the formula: 


Ediss[R{t)] 



dt'j[R{t')]R{t'f, 


(17) 


until time t corresponding to the relative distance R{t) in 
the entrance channel. Figure [T51 shows the accumulated 
dissipation energy as a function of R for the ®®Zr -|- ^^^Sn 
system at the same Ec.m. as in Fig. [TH A zoom around 
A = 13 fm is inserted in the figure to focus on a low Adiss 
region. As expected from the property of the extracted 
friction coefficient, the dissipation energies are zero at 
large A and then monotonically increase as the colliding 
nuclei approach and overlap each other. Irrespective of 
humps and peaks of the extracted friction coefficients, 
the dissipated energies show a smoothed behavior. Con¬ 
cerning the dependence of dissipation energy on Ec.m. of 
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FIG. 17. (Color online) Schematic picture to define the quan¬ 
tities used in the analysis: Frozen density potential bar¬ 
rier Vfd, fusion threshold energy Sthres, extra-push energy 
l^extra = Ethrea — VpD, and increase in potential AV for the 
®®Zr -I- ^^'’Sn system. 

TDHF trajectory, the dissipation energy becomes larger 
at larger Ec.m. because of larger relative velocity entering 

in the integral in Eq. ( 113 - 

In the following, we will use the dissipation energy 
Ediss accumulated until i?min as the physical quantity 
to analyze the fusion hindrance in the entrance channel 
of fusion reactions. 

D. Origin of extra-push energy 

In Sec. IIII A[ we have observed the fusion hindrance 
phenomenon by TDHF simulations for the entrance chan¬ 
nel of central collisions in heavy systems. We have shown 
in Fig. 13] the extra-push energies estimated from TDHF 
simulations for several heavy systems and reasonably re¬ 
produced the extra-push energies deduced from experi¬ 
mental observations of evaporation-residue cross sections. 
In this subsection, using the properties of the extracted 
potentials and friction coefficients in heavy systems pre¬ 
sented in the previous subsections, we will analyze the 
origin of the fusion hindrance and appearance of the 
extra-push energy. 

Before showing the result of the analysis, we first sum¬ 
marize the quantities in a schematic picture in Fig. [T7| 
that will be used in the analysis. They are the fu¬ 
sion threshold energy Ethres and the frozen density po¬ 
tential barrier Vfd, indicated by the arrows, for the 
®®Zr-|-^^^Sn system. We have defined the extra-push 
energy by TDHF as the difference between the fusion 
threshold energy and the frozen density potential bar¬ 
rier, Eextra = A'thres — VpD • Here, we introduce a new 
quantity, increase in potential energy AH, that is the dif¬ 
ference between the extracted potential at Ethres at Rmin 
and the frozen density potential barrier. Since the ex¬ 
tracted potential increases monotonically as R decreases 
in the heavy systems, AV becomes maximum at 
We will also use in the analysis the accumulated dissipa¬ 
tion energy E^iss defined in Eq. (EZD. This accumulated 


dissipation energy is computed from the friction coeffi¬ 
cient extracted at Ethres- As shown in Eig.[Tni E^iss also 
monotonically increases as R decreases. 

Erom the energy conservation, we have at i? > i?niin 

Ec.m. = Aki„(A) + H(i?) + Ediss(A), (18) 

where Ekin = H^/2/i(i?) is the collective kinetic energy 
and all the energies in the right-hand side can be com¬ 
puted by TDHF and by the DD-TDHF method and we 
confirm that the energy conservation holds. Subtracting 
Vfd from the both sides of Eq. (fTSl) and using the quan¬ 
tities Eextra = Athres “ V^D and AV = V^(i?min) “ HpD, 
we derive the following relation at Ec.m. = Athres at 

R — .^min- 

Aextra — E]cm{Rmin )+AV{Rmm) + ^diss (^min ) (19) 

As shown in Eig. jS] the absolute value of the momen¬ 
tum |P| is expected to monotonically decrease and to be 
smallest at i? = Rmin in heavy systems, so is Ekin- We 
have checked that Ekin at i? = i?min in all the systems is 
less than 0.6 MeV and small enough, compared with AV 
and Ediss- Therefore, we approximate the extra-push en¬ 
ergy as the sum of AH and Ediss, 

Aextra « AV{Rmm) + -^diss (-^min ). (20) 

That is, we can identify the origin of the extra-push en¬ 
ergy in the entrance channel as AH and/or Ediss ex¬ 
tracted from the DD-TDHF method. 

Figure [TH| shows each contribution to the extra-push 
energy for the same systems as in Fig. [31 The red 
filled circles denote the contribution of the increase in 
potential AH to the extra-push energy, while the blue 
filled triangles denote the contribution of the accumu¬ 
lated dissipation energy Ediss- In the 90Zr-|- 
and ^°°Mo-V ®^Mo systems, AH is smaller than Ediss- In 
the rest of the systems considered here, AH is greater 
than Ediss- As the charge product of the system in¬ 
creases, the ratio of AH to Ediss increases. From this 
figure, we can conclude that, as the system becomes heav¬ 
ier, the contribution of the increase in potential energy to 
the extra-push energy in the entrance channel is more im¬ 
portant than that of the accumulated dissipation energy. 
Dissipation plays an important role, but is not sufficient 
to explain the amount of the extra-push energy in our 
analysis with the DD-TDHF method. 


IV. CONCLUSIONS 

The synthesis of superheavy elements in laboratories 
is extremely challenging because of an extremely low 
production rate. One of the reason for a low produc¬ 
tion late is the presence of the fusion hindrance. This 
gives strong competition between compound nucleus for¬ 
mation and quasifission. In this article we address this 
point by the selfconsistent microscopic reaction theory. 
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FIG. 18. (Color online) Contributions of the increase in potential AV and of the accumulated dissipation energy i?diss to the 
extra-push energy for the same systems as those in Fig. [S] 


TDHF. We have shown the applicability of TDHF dy¬ 
namics for fusion reactions involving heavy nuclei by re¬ 
producing the extra-push energies estimated from ex¬ 
perimental evaporation-residue cross sections in heavy 
systems. Then, we have analyzed the property of the 
nucleus-nucleus potentials for the entrance channel of 
fusion reactions extracted from the DD-TDHF method, 
which combines TDHF dynamics with a classical equa¬ 
tion of motion. In heavy systems, the extracted nucleus- 
nucleus potentials show the following properties: (i) Or¬ 
dinary Coulomb barrier in the extracted potential disap¬ 
pears and the potential monotonically increases as the 
relative distance decreases, (ii) The dependence of po¬ 
tential on the energy used in the TDHF simulations is 
not pronounced. These properties are not seen in the 
potentials in light- and medium-mass systems and also 
in the DC-TDHF method. We have shown that the in¬ 
crease in potential is a reflection of the appearance of the 
conditional saddle inside the contact configuration of col¬ 
liding nuclei on an adiabatic potential surface. By ana¬ 
lyzing the isotope dependence of the extracted potentials, 
we have found that, as nucleon transfers towards charge 
equilibration during the collision becomes small, increase 
in the extracted potential from the frozen density poten¬ 
tial becomes large. Extracted friction coefficients how¬ 
ever show energy dependence, and have a sharp peak 
near the relative distance at which the extraction stops. 
This indicates strong reorganization of the internal struc¬ 
ture near that distance. Using these properties of the ex¬ 
tracted nucleus-nucleus potential and energy dissipation, 
we have analyzed the origin of the extra-push energy in 
heavy systems. The analysis has shown that more con¬ 
tribution comes from the increase in potential energy to 
extra-push energy than that from the accumulated dis¬ 
sipation energy as the size of the system is heavier. We 
have concluded that the increase in potential is more im¬ 
portant than dissipation for understanding the origin of 
the fusion hindrance. 

Although the present study is limited to central colli¬ 
sions, which enable us to perform the DD-TDHF method 
with one-dimensional equation of motion, we have shown 


the feasibility of the present analysis for the fusion hin¬ 
drance. Possible extensions to non-central collisions in¬ 
volving finite angular momenta will be interesting to un¬ 
derstand the role of angular momentum for the fusion 
hindrance. Furthermore, the orientation dependence of 
the dynamics due to large nuclear deformation will be 
important for systems with, for example, actinide nuclei, 
which have been and will be used in the so-called hot 
fusion reactions in the synthesis of superheavy elements. 
The fission process, as an inverted process of fusion, is 
one of the future subject to apply the present method to 
extract information on energy dissipation. 
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Appendix A: Criterion of determining the neck 
position 

To construct two-body kinematics in the DD-TDHE 
method in Sec. KSl we need to separate the total sys¬ 
tem into two subsystems at the neck position. In the case 
of central collisions along x-axis, the neck position can be 
specified hy x = xq. Previously, we have determined the 
neck position where the densities of the projectile and 
target cross and checked that this definition is good for 
light- and medium-mass systems at any energies [31^ , 
hereafter called criterion (i). However, we realize that 
this criterion is inappropriate for some cases, especially at 
small relative distances in heavy systems. Then, we test 
the following criteria to better determine the neck posi¬ 
tion: (ii) The position where the total one-dimensional 
density p{x,y = yoiO) becomes minimum at yo = Ofm. 
(iii) The same as criterion (ii) but at yo = 2 fm. By look¬ 
ing at the two-dimensional density profiles, we find that 
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the criterion (iii) is most reasonable and stable for the 
cases at low energies in most heavy systems considered 
in this article. Therefore, we use the criterion (iii) to 
determine the neck position at low energies. For higher 
energies, we find that the criterion (i) is better. 


Appendix B: Ambiguity of determining the fusion 
threshold energy 

We have explained how to determine the fusion thresh¬ 
old energy in Sec. Ill PI When we determine the fusion 
threshold energy, we have encountered an ambiguity to 
define a compact shape in some systems. In the follow¬ 
ing, we will show that the conclusion of our study, i.e., 
the increase in potential is more important than the accu¬ 
mulated dissipation energy for understanding the origin 
of the fusion hindrance, is not be affected by this am¬ 


biguity. To do so, we take the ®°Zr-|-^^^Sn system as 
an example. In this system, we indeed define the fusion 
threshold energy as i?c.m. = 234.5 MeV. The reason is 
that Q 20 at Ec.-m. = 234 MeV is 15% larger than <520 at 
higher energies (see Fig. [1]) and we regard the configura¬ 
tion at i?c.m. = 234 MeV as not a compact shape. There¬ 
fore, we exclude the case at Ec.m. = 234 MeV from the 
fusion events. However, one might consider this case as 
still a compact shape, if one compares this with the case 
at Ec.m. = 233.5 MeV, in which the system reseparates. 
In order to check how this ambiguity affects the result, 
we change the fusion threshold energy from 234.5 MeV 
to 234 MeV in this system and repeat the same calcula¬ 
tion. We find that the increase in potential and accu¬ 
mulated dissipation energies are changed from 8.3 and 
7.9 MeV for Flthres = 234.5 MeV to 8.0 and 7.7 MeV for 
i^thres = 234 MeV, respectively. The ratio of AV/Ediss 
does not change significantly with the change of Flthres 
and we therefore hold the conclusion as it is. 
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